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Non-trivial geometrical effects in relativistic central collisions of deformed nuclei 
are studied using a simple version of optical Glauber model. For very small impact 
t — ' parameters large centrality and eccentricity fluctuations are observed. In very high 
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multiplicity collisions of oblate nuclei (e.g. Au 197 and Cu) a significant fraction of 
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events with elliptic flow strength V2 dependent on oblateness B 2 is predicted. 
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I. INTRODUCTION. 

In relativistic heavy-ion collisions a simple relation between centrality of collisions c(N) 
and average participant accent of interring nucleons < <p > is usually expected to 
be valid. According to rather general assumptions [1[ centrality c(N) is directly related to 
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impact parameter b as c(N) ~ 7tb 2 /ai ne i. At the same time the participant eccentricity 
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ep(6) = E(^-^)/E(^+l/?) C 1 ) 



(where Xi and are positions of interacting nucleons in transversal plane) is also singly 
dependent (within fluctuations) on the size of impact parameter b. Thus measured elliptic 
flow strength v 2 is expected to rise from very small values in most central collisions 

proportionaly to increasing values of average eccentricity v% ~ k- < e p (b) >. 

However these simple relations are not valid for collisions of deformed nuclei. The orien- 
tation of a deformed nucleus (e.g. relative to beam axis) has a direct influence on centrality 
and eccentricity of collisions at a given fixed value of impact parameter b. This is slightly 
worrying since the elliptic flow strength measured in collisions of nuclei having non-zero 
deformation /? 2 (e.g. Au, Cu) has been interpreted Q] assuming these nuclei to be spherical. 

In the following sections we present a preliminary study of non-trivial variations of initial 
eccentricity at given collision centralities in relativistic collisions of deformed nuclei. 
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Fig.l: Predicted shapes of Pb 207 , As 75 (/3 2 = -0.25) and Ho 165 {(5 2 = 0.3) nuclei [9]. 
II. GEOMETRICAL EFFECTS IN COLLISIONS OF DEFORMED NUCLEI 

Although effects of nuclear deformation have been carefully studied in fusion reactions at 
coulomb barrier [5J the role of nuclear deformation in relativistic heavy ion collisions have 
been considered scarcely so far. It had been pointed out 6i] that highest energy densities of 
QCD matter could be created in central collisions of longitudinaly oriented heavy deformed 
nuclei and elliptic flow in such collisions had been discussed. However since the heaviest 
nuclei suitable for relativistic collisions experiments are prolate the elliptic flow in central 
collisions of oblate nuclei have escaped the attention of heavy ion community so far. 

In this contribution we point out that oblate ground-state nuclear deformation of stable 
nuclei can have a significant influence on measured elliptic flow values in central collisions 
since it directly influences fluctuations and average values of initial eccentricity e p at given 
collision centralities. 

Let us consider angle 9 between the beam direction and the main axis of a deformed 
nucleus (e.g. Ho-165). For 9 ~ nucleus is longitudilay polarized and for 9 ~ tt/2 we have 
transversal polarization. Even for very small impact parameters b ~ (central collisions 
geometrically) the orientation of a deformed nucleus relative to beam axis influences directly 
participant eccentricity 

T* x 2 - v 2 
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and thus final strength of the elliptic flow observed. Also centrality of collision depends on 
the orientation of a deformed nucleus at given impact parameter. This happens not only 
due to changes in the number of participating nucleons (a simple overlap effect) but also due 
to that fraction of secondary particles multiplicity which is proportional to total number of 
binary nucleon-nucleon collisions Ml. 

n 

Thus geometric relation between centrality and impact parameter [1[ is not valid for 
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collisions of deformed nuclei and initial eccentricity e p is not simply related to collision 
centrality as it is usually assumed. This can in principle modify interpretations drawn from 
the eliptic flow strength t> 2 at given centrality c(N) in Au+Au collisions at RHIC. 



III. A SIMPLE OPTICAL GLAUBER SIMULATION 



In order to test quantitatively the effects of nuclear deformation on centrality and ec- 
centricity in heavy ion collisions a simple Optical Galuber Model (OGM) calculation for 
deformed nuclei had been performed. Here is a very short description of the simulation: 

In the first step a deformed Woods-Saxon [s| density 
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with deformation parameters (3%, 04 taken from 9| for various nuclei (Au,Cu,H°>Pb,Si,Ca) 
had been used to generate transversal projections p T (x,y) = J p(r)dz of nucleon density. 
Binning of p T (xi, yj) histograms was 0.1/m and longitudinal integration of deformed Woods- 
Saxon density p(r) was performed in 0.05/m steps. For deformed nuclei separate projections 
p T [9, 4>] had been generated at fixed values of angle 6 and azimuthal angle </>. 

| Pb-20B <b2^0.0) | | Cu-63 (b2=Q.16) | 




Fig. 2: Projected nucleon densities p T (x,y) for Pb 208 , ^Cu 63 , ^Cu 65 , ^Au 197 (/? 2 from 



In order to simulate also non-central collisions transversal projections p T (x,y) had been 
shifted for selected values of impact parameter b = (b x , b y ) and total number of nucleon- 
nucleon collisions N co n had been evaluated as 

N coU = Ncouixi, Vj ) = £>aw ' Pl{xi - yj - b y /2) * pl{ Xi + 6 x /2, Vj + b y /2). (4) 

The cross section g™n was fi xe d at 42mb. 

Using the geometrical information stored in histograms pj(xi,yi) and p%{xi,yj) transver- 
sal projected baryon density pb(x, y) and number of participants N part were simply evaluated 
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and found to be in good agreement 1 with results of Monte-Carlo Glauber model simulation 
[lo| as shown in Figs. 9-10 (Appendix). 

Eccentricity €nn of the interacting volume had been evaluated as 

~ Y,i,jN coll {x i ,y j ) ■ [x\-y]) 

' VV >:,,A,,//(.'',// ; )-(.rr • ///) 1 j 

in CMS coordinates of the overlaping zone. Numerical values of eccentricity e^N were found 

to be in reasonable agreement with Au+Au eccentricities evaluated in [3J (see Fig. 11 in 

Appendix). Selected results of these calculations for Au,Cu,H° ; Pb,Si nuclei are presented 

in the following sections. 



IV. COLLISIONS OF SPHERICAL + DEFORMED NUCLEI 

For the sake of simplicity let us start with polarized collisions of a deformed nucleus e.g. 
a prolate Ho-165 (fa ~ 0.3) with spherical projectile Pb-207. 
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Fig. 3: Eccentricity [n-n collisions] plot for Pb+Pb and polarized ^Ho+Pb systems. 

If the spin^ 7 / 2 ) of Ho-165 nucleus is orthogonal to beam axis (T-polarized) the participant 
eccentricity e p will be non-zero even for b = its value roughly corresponding to non-central 
Pb+Pb collisions at b = 4 fm (see Fig. 3). Increasing impact parameter values (b x ,b y ) in 
directions parallel (b y ) and orthogonal (6^) to Ho 165 spin have different influence on the 
eccentricity. For increasing (b x , 0) the eccentricity further increases while for increasing 



1 Spherical shape for Au-197 with R = 6.38fm and diffusivity a = 0.53fm was used for this comparison. 
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(0, by) the eccentricity first decreases to zero at b y ~ 4 fm and then it increases as it can be 
seen from Fig. 3 (data points are evaluated for impact parameters increasing in steps lfm). 

For longitudinal polarization of holmium nuclei (relative to beam direction) the initial 
eccentricity tends to zero in the most central collisions. This happens because eccentricity 
fluctuations |TjJ are not taken into account in our simple optical Glauber calculation. It is 
also visible from Fig.3 that central collisions of L-polarized Ho-165 with Pb have significantly 
higher multiplicity compared to central T-polarized Ho+Pb collisions (with 6 = 0). Thus 
the highest multiplicity collisions of prolate nuclei correspond to longitudinal orientations of 
deformed nucleus due to the highest number of nucleon-nucleon interactions in this case. 

This means that in order to study the elliptic flow in central collisions of transversaly ori- 
ented prolate nuclei with spherical projectiles one needs to select events with lower measured 
multiplicity then maximal. In such sample of events non-central collisions of longitudinaly 
oriented prolate nuclei will get mixed together with central collisions (b ~ 0) of transversaly 
oriented prolate nuclei of similar multiplicity (see Fig.3). One can try to distinguish these 
two types of events via measurement of spectator energy in ZDC calorimeters and to perform 
target (or beam) polarization-dependent experiments (see [3] for polarized Ho-165 target). 



| Au+Pb Eccentricity vs. Centralily | Ca40+Si29 Eccentricity vs. Centrality __) 




Ncoll. [ n-n collisions ] Ncoll. [ n-n collisions ] 

Fig. 4: Eccentricity plots for polarized ^Au 197 +Pb 207 and ^Si 29 +Ca 40 collisions. 

Fortunately the existence of nuclei with oblate (/3 2 < 0) deformation (Au,I,Ga,As,Si,Al) 
provides us with another possibility in this direction. For oblate nuclei very high multiplicity 
collisions correspond to transversal polarization of nuclear spin relative to the beam direction 
due to higher number of nucleon-nucleon collisions in this case. Thus significant fraction 
of very central collisions of oblate nuclei should exhibit non-zero elliptic flow t>2 values 
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proportional to deformation parameter /3 2 < 0. In Fig. 4 we show results of our calculations 
for Au 197 +Pb 207 and Si 29 +Ca 40 collisions. A small difference in the number 2 of nucleon- 
nucleon interactions in the most central collisions with longitudial and transversal spin 
polarizations is obvious. In the next section we show that the influence of nuclear ground- 
state deformation on the eccentricity fluctuations in Au+Au and Cu+Cu collisions (studied 
at RHIC) is expected to be more significant. 

V. COLLISIONS OF TWO DEFORMED NUCLEI 

In collisions of two deformed nuclei various geometrical configurations are possible. Let 
us denote Q\ and 82 to be angles of deformed nuclei main axes relative to beam direction 
and 0i and 02 to be main axes azimuthal angles in LAB. For unpolarized beams (targets) 
azimuthal angles are independent and nuclei collide at random orientations of 0i,02- A 
trivial consideration shows that probability of having a collision with deformed nuclei axes 
oriented at relative azimuthal angle < A0 < 7r/2 is constant: P(A0) = const. 




Polar angles 8\, 82 are also mutually independent but not random and probability P(8\, 62) 
of a collision of two nuclei oriented at angles < Q\ < it and < 9 2 < n is 

P(9 1 ,9 2 ) =sin(0i)sin(0 2 )/4 (6) 

which says that longitudinal-longitudinal oriented (Long-Long) collisions of deformed nuclei 
are rather rare while transversal (Trans- Trans) collisions {9\ ~ 82 ~ vr/2) are more frequent. 

One should keep in mind that this applies to collisions at any value of geometrical impact 
parameter occuring with probability P{b < b max ) = c ■ b 2 max . Therefore central {b rs 0) 



2 One should always keep in mind the existence of fluctuations in the number of n-n collisions. 
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Long-Long collisions are very rare! Nevetherless, just such collisions of prolate nuclei can be 
selected |6| from the sample of events according to the centrality (high-multiplicity) criterion 
and analyzed. 
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Fig. 6: Densities N co u(x, y)/fm 2 for Au^+Au^ at 6 x =0fm/3fm and Eccentricity [Ncoll] plot. 



As can be seen from Fig. 6 Trans- Trans collisions of oblate nuclei with parallel orientation 
of nuclear spins (A<f) = 0) are the highest multiplicity collisions for nuclei with oblate 
geometry. They can be localized in the sample of events based on the collision centrality 
and elliptic flow v 2 which is expected to be increasing with oblateness strength -^ 2 - 
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Fig.7: EccfNcoll] plot for Cu^+Cu 63 and Cu 65 +Cu 65 assuming (3 2 = 0.16 and -0.15 9] 



If stable copper isotopes are deformed as predicted in [9] then most central collisions of 
Cu 63 should exhibit significantly different elliptic flow strength t> 2 compared to most central 
collisions of Cu 65 +Cu 65 (predicted ft 2 = -0.15). This means that elliptic flow measured in 
most central collisions of a given isotope can be used to determine f3 2 sign of its ground-state 
deformation. This might hold also for spin-zero nuclei e.g. Sig+ (predicted (3 2 = —0.48 
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Situation is however slightly more complex than indicated in Figs. 6 and 7. In the case 
of Trans- Trans collisions (9i ~ 9 2 ~ 7r/2) of deformed nuclei there is a freedom in relative 
azimuthal orientation A0 of nuclear spins. For A<f) m n/2 (orthogonal spins) elliptic flow 
should tend zero and for A0 » (parallel spins) participant eccentricity and thus also the 
observed elliptic flow strength v 2 should be maximal (whatever the sign of j3 2 is). 

| Si29+Si29 Eccentricity vs. Centrality | 
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Fig.8: Eccentricity [Ncoll] plot for T Si 29 + T Si 29 assuming /3 2 =-0.3 and /3 4 =0.13 



a. 



This is clearly visible from our calculations in Figure 8, where dashed arrow indicates 
increase of eccentricity for transversally polarized 9\ = 9 2 = 7r/2 collisions of Si 29 +Si 29 when 
azimuthal angle difference of spins changes from A0 = 7r/2 (denoted as Trans- Trans RT) 
to parallel spins orientation A<p = 0. Therefore in geometrically most central collisions 
(zero impact parameters b) the participant eccentricity depends on polar angles Q\ , 9 2 and 
on relative azimuthal angle difference A0: 



e p [bnO] = e p (9 1 ,9 2 ,A<j>) (7) 

For random relative azimuthal angles A0 a sample of very high multiplicity events will 
contain mixture of collisions with different A</> and various eccentricities. This applies also 
to Au+Au collisions studied at RHIC. Collisions of longitudinaly+transversaly polarized 
Si 29 have lower maximal number of binary collisions N co u w 80 (at 6=0 fm) with eccentricity 
£nn ~ 0.15 (see Fig.8). Data points in all our Eccentricity[Ncoll.] plots correspond to 
impact parameters h x and b y increasing in steps 0,1,2,3,4,5,6 fm. Calculated number of 
binary collsions decreases with increasing impact parameters (b x ,0) and (Q,b y ). 
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As a matter of completeness we mention here also the possibility to collide prolate nuclei 
(fa > 0) with oblate nuclei (fo < 0) such as Ho 165 +As 75 or Cf 251 +Au 197 . In this case 
collisions with highest multiplicities will be Long- Trans polarized and they should exhibit 
elliptic flow strength < v 2 > rising with (5 2 deformation of the prolate nucleus. 

VI. DISCUSSION OF RESULTS 

It has been shown that in collisions of transversaly polarized prolate nuclei (e.g. Ho-165) 
with spherical projectiles (e.g. Pb-207) a significat elliptic flow will be generated even for 
zero impact parameters. Such collisions however will have similar multiplicity (centrality) 
as non-central collisions with longitudinaly polarized prolate nuclei. Careful polarization- 
dependent studies and precise measurement of spectator energy can distinguish different 
orientations of prolate nuclei relative to beam axis in such collisions. Highest-multiplicity 
collisions of prolate nuclei (such as Ho-165) correspond to longitudial orientations of prolate 
nuclei. Elliptic flow strength v 2 will thus tend to zero (within fluctuations) in such collisions. 

For oblate nuclei the situation is different. Highest multiplicity collisions of oblate nucleus 
with a spherical projectile correspond to transversal polarization of oblate nucleus relative 
to beam axis. Additionally if both nuclei colliding are prolate highest multiplicity collisions 
will tend to be collisions of transversally polarized nuclei with parallel azimuthal orientation 
of nuclear deformation axes. Average value of elliptic flow strength v 2 measured in very high 
multiplicity (VHM) collisions of oblate nuclei will thus contain two contributions: 

< v 2 > VHM =< v 2 > fluct + <v 2 > P2 (8) 

The first contribution < v 2 > fi ucL corresponds to eccentricity fluctuations which are present 
also in the highest multiplicity collisions of spherical nuclei. The second contribution is 
proportional to geometrical deformation of oblate {(3 2 < 0) nucleus. Quantities < v 2 > fi uct 
and < v 2 >(i 2 can in principle be disentangled if elliptic flow for highest multiplicity collisions 
of spherical nuclei with similar multiplicity is measured < v 2 > l/HM =< V2 >f luct . This 
approach can be attempted for two isotopes of a given element (e.g. Si 29 and Si 30 ) or for 
two nuclei (spherical+oblate) with similar number of nucleons e.g. Pb and Au (see Fig. 6). 
One should keep in mind still that < v 2 >p 2 contains averaged values of t> 2 (A0) due to 
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almost 3 random relative azimuthal orientations A0 of nuclear spins in such highest multi- 
plicity events (for unpolarized beams and targets). 



VII. CONCLUSIONS 

A simple simulation based on optical limit approximation of the Glauber model [13] shows 
that very high multiplicity collisions of oblate nuclei at high energies should exhibit non-zero 
elliptic flow strength dependent on fa deformation parameter of the oblate nucleus. This 
provides us with the possibility to study the elliptic flow phenomenon in the most central 
ultra-relativistic collisions of oblate nuclei. The equation of state of QCD matter can thus 
be investigated in the highest multiplicity collisions. The energy dependence [ijj] of elliptic 
flow strength v 2 in these very high multiplicity (VHM) central collisions of oblate nuclei can 
reveal changes in the equation state of hadronic matter (e.g. kink-like signatures [if] ]) 
due to the expected phase transition of QCD matter. 

For a precise investigation to what extent the oblateness of Au 197 nucleus influences 

n 

interpretations drawn from elliptic flow measurements at RHIC |4[ a full MC simulation 
based on Glauber model taking into account ground-state deformation of Au 197 nucleus 
together with all possible orientations 0i,9z,<t>i, 4>2 of colliding nuclei is probably necessary. 
A comparison of experimental results from Au+Au collisions with collisions of spherical 
nuclei (e.g. Pb 207 ) at the same energy range and the study of Cu+Cu collisions with the 
second stable copper isotope at RHIC could be very useful. Elliptic flow studies in collisions 
of deformed nuclei at LHC would provide us with new and possibly exciting results. 



VIII. APPENDIX: COMPARISON WITH MC GLAUBER CALCULATIONS. 

For a comparison with results of full Monte-Carlo Glauber calculation [lCj we present 
here impact parameter dependencies of some quantities obtained with our simple optical 
Glauber model (OGM) assuming Au 197 nucleus to be spherical (R = 6.38fm; a = 0.53). 

The impact parameter dependence of the number of binary nucleon-nucleon collisions 
together with number of participants and spectators calculated in our siplified OGM model 

3 The influence of Acj) on collision centrality at b « is rather small (see Fig. 8 for Si+Si collisions). 
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Comparison with other MC \ 
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Fig. 9: Centrality dependence of binary nucleon-nucleon collisions and participants in Au+Au. 

is shown in Figure 9. Data points from full MC Glauber simulation [h]] are shown for a 
comparison. 

In Figure 10 we show the maximal transverse baryon density in the overlapping zone and 

maximal number of nucleon-nucleon collisions per fm 2 in transversal plane for spherical Au 

collisions. Our simple OGM calculations are in reasonable agreement with results of full 

□ 

MC Glauber simulation 10 . 
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Fig. 10: Maximal baryon density x [&] and maximal collision density p^coU ft] i n Au+Au. 

Impact parameter dependence of our calculated (OGM) eccentricity e^v for Au+Au 
collisions (assuming spherical Au 197 nucleus) together with eccentricity values calculated by 
R.S.Bhalerao et.al. [3( is shown in Figure 11. 
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Comparison of Eccentricity Au+Au | 
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Fig. 11: Impact parameter dependence of (OGM) eccentricity £nn in comparison with results 
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